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I .  INTRODUCTION 


Simpson's  One-Third  Rule  in  composite  form  is  one  of  the  most 
frequently  used  of  the  simpler  quadrature  formulas  for  equidistant 
ordinates.  The  analysis  of  error  in  the  rule  has  been  extensive,  and 
relations  are  available  for  estimating  the  error.1  The  rule  is 
applicable  only  for  an  even  number  of  equidistant  subdivisions  of  the 
range  of  integration,  and  the  coefficients  of  the  ordinate  values  away 
from  the  endpoints  differ  by  a  factor  of  two. 

The  purpose  of  this  report  is  to  present  a  composite  formula 
generated  from  Simpson's  Rule  that  applies  for  either  an  odd  or  even 
number  of  subdivisions  of  the  range  of  integration,  and  whose 
coefficients  of  the  ordinates  are  unity  except  near  the  ends  of  the 
integration  range.  Having  coefficients  of  unity  over  most  of  the 
integration  range  is  an  advantage  when  the  ordinates  correspond  to 
measurements  with  errors,  because  the  errors  are  weighted  equally  in  the 
integration  process.  Amplification  of  roundoff  noise  in  computed 
ordinates  will  be  minimized  also.  However,  to  obtain  these  advantages, 
two  additional  ordinates  must  be  obtained  and  utilized. 


II.  DEVELOPMENT  OF  THE  FORMULA 
Simpson's  One-Third  Rule  in  composite  form  is: 

I  -  (h/3)(y0  +  Ayi  +  2y2  +  ....  +  2yn_2  +  4yn_1  +  yn)  +  E 
where:  I  =  value  of  integral 

y  ,  ..y„  =  ordinate  values  at  intervals  of  h  on  the  abscissa 
E  =  error  term 

Consider  a  range  of  integration  divided  into  an  odd  number  of 
equidistant  intervals  of  width  h.  The  rule  can  be  applied  to  the  range 
of  an  even  number  of  intervals  with  ordinates  of  yQ  through  yn_^  •  The 

rule  also  can  be  applied  to  the  range  with  ordinates  of  v-^  through  y  . 

If  the  two  results  are  added  the  sum  is  twice  the  value  of  the  desired 
integral  less  the  values  of  the  integral  for  the  intervals  between  yo  and 

y^  and  between  yn_j  and  yn .  If  a  value  for  the  ordinate  at  the  midpoint 

of  each  of  these  end  intervals  can  be  obtained,  then  Simpson's  Rule  can 
be  applied  to  obtain  the  integral  values  for  the  end  intervals.  Adding 
these  values  to  the  previously  obtained  sum  from  the  application  of  the 


J.  B.  Scarborough,  Numerical  Mathematical  Analysis,  4th  Edition,  The 
Johns  Hopkins  Press,  Baltimore,  Maryland,  1958. 

R.  W.  Hamming,  Numerical  Methods  for  Scientists  and  Engineers,  McGraw- 
Hill  Book  Company,  New  York,  New  York,  1962. 
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composite  formula,  the  following  result  for  the  full  range  of  integration 
for  an  odd  number  of  subdivisions  of  width  h  is  obtained: 

21  -  (h/3) (y0  +  4yx  +  2y2  +  4y3  + - +  2yn_3  +  4yn_2  +  yn-1) 

+  (h/3) (y  x  +  4y2  +  2y3  +  ...  +  4yn_3  +  2yn_2  +  4yn_1  +  yn) 

+  (h/6) (y0  +  4y1/2  +  yx)  +  (h/6)  (yn_1  +  4yn-i/2  +  yn>  +  E 

After  combining  terms,  the  formula  becomes: 
i  =  h  £ (1/ 4)y0  +  (i/3)y1/2  +  (n/i2)yi  +  y2  +  y3  +  ••••  +  yn_3  +  yn_2 

+  ( 1 1/1 2)yn_1  +  (l/3)yn_1/2  +  (l/4)yn  ]  +  E  (1) 

All  coefficients  except  for  the  end  points  and  the  midpoints  of  the 
end  intervals  are  11/12  or  one. 

Now  consider  a  range  of  integration  divided  into  an  even  number  of 
equidistant  intervals  of  width  h.  Here  Simpson's  Rule  can  be  applied  to 
the  entire  range  of  integration.  It  can  also  be  applied  to  the  range 
beginning  with  the  ordinate  y^  and  ending  with  the  ordinate  yn_^>  If  the 

results  are  added  the  sum  as  before  is  twice  the  desired  value  for  the 
integral  less  the  integral  values  for  the  intervals  from  yQ  to  y3  and 

from  yn_i  to  yn .  Assuming  that  the  ordinate  for  the  midpoint  of  each  end 

interval  is  available,  the  integral  values  for  the  end  intervals  can  be 
obtained  by  applying  Simpson's  Rule,  and  the  formula  for  the  entire 
integration  range  becomes  as  follows: 

•  21  =  (h/3) (y0  +  4Yl  +  2y2  +  ....  +  2yn_2  +  4yn_1  +  yn) 

+  (h/3) (y x  +  4y2  +  2y3  +  -  +  2yn_3  +  4Yn_2  +  yn-1) 

+  (h/6)(y0  +  4Yi/2  +  yx)  +  (h/6)  (y^  +  4yn_1/2  +  yn)  +  E 

After  terms  are  combined  the  formula  for  I  is  found  to  be  identical 
to  that  presented  in  Equation  1 .  Therefore  Equation  1  applies  for  both 
even  and  odd  n . 

Because  the  integral  for  the  two  end  intervals  is  more  accurately 
determined  through  use  of  one-half  width  intervals  in  the  computation, 
the  error  for  the  full  integration  range  will  be  somewhat  less  than  that 
calculated  for  Simpson's  Rule  for  n  subdivisions.  However,  the  reduction 
is  small,  and  becomes  less  as  the  number  of  subdivisions  is  increased, 
since  the  proportion  of  the  range  determined  with  less  error  becomes 
smaller  as  the  number  of  subdivisions  increases.  Therefore  error 
estimates  for  the  formula  in  Equation  1  are  made  as  for  Simpson's 
Composite  Rule,  ignoring  the  one-half  interval  ordinates. 
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In  summary,  the  quadrature  formula  presented  in  Equation  1  applies 
for  an  odd  or  even  number  of  equidistant  subdivisions  of  the  range  of 
integration,  and  provides  coefficients  of  unity  or  near  unity  except  fo 
the  two  ordinates  at  each  end  of  the  integration  range.  The  error 
estimate  can  be  made  as  is  done  for  Simpson's  Composite  Rule.  However, 
to  apply  the  formula,  two  additional  ordinates  must  be  available:  fhos 
at  the  midpoint  of  the  intervals  at  each  end  of  the  integration  range . 
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1  Agbabian  Associates 
ATTN:  M.  Agbabian 

? 50  North  Nash  Street 
El  Segundo,  CA  90245 

1  The  BDM  Corporation 

ATTN:  Richard  Hensley 

P.0.  Box  9274 
Albuquerque  International 
Albuquerque,  NM  87119 

1  The  Boeing  Company 

ATTN:  Aerospace  Library 

P.0.  Box  3707 
Seatle,  WA  98124 

1  Coodyear  Aerospace  Corp 

ATTN:  R.  M.  Brown,  Bldg  1 

Shelter  Engineering 
Litchfield  Park,  AZ  85340 

5  Kaman  AviDyne 

ATTN:  Dr.  N.P.  Hobbs  (4  cys) 

Mr.  S.  Criscione 
83  Second  Avenue 
Northwest  Industrial  Park 
Burlington,  MA  01830 

3  Kaman  Sciences  Corporation 
ATTN:  Library 

P.  A.  Ellis 
F.  H.  Shelton 

1500  Garden  of  the  Gods  Road 
Colorado  Springs,  CO  80907 


No  .  of 

Co  p i e  s  Organizat Ion 

1  Kaman  Sciences  Corporation 
ATTN:  Don  Sachs 

Suite  703 

2001  Jefferson  Davis  Highway 
Arlington,  VA  22202 

1  Kaman-TEMPO 
ATTN:  DAS I AC 

P.0.  Drawer  0Q 
Santa  Barbara,  CA  93102 

1  Kaman-TEMPO 

ATTN:  E.  Bryant,  Suite  l'L-1 

715  Shamrock  Road 
Bel  Air,  MD  21014 

1  Lockheed  Missiles  &  Space  Co. 

ATTN:  J.  J.  Murphy,  Dept.  81-11 

Bldg.  154 
P.0.  Box  504 
Sunnyvale,  CA  94086 

1  Martin  Marietta  Aerospace 
Orlando  Division 

ATTN:  G.  Fotieo 

P.0.  Box  5837 
Orlando,  FL  32805 

2  McDonnell  Douglas  Astronautics 

Corporation 

ATTN:  Robert  W.  Halprin 

Dr.  P.  Lewis 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 

2  The  Mitre  Corporation 
ATTN:  Library 

J.  Calligeros,  Mail 
Stop  B-150 
P.0.  Box  208 
Bedford,  MA  01730 

1  Pacific  Sierra  Research  Corp 
ATTN:  Dr.  Harold  Brode 

1456  Cloverfield  Boulevard 
Santa  Monica,  CA  90404 
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DISTRIBUTION  LIST 


No .  of 

Copies  Organization 

2  Physics  International  Corp 
ATTN:  Technical  Library 

Mr  .  Fred  Sauer 
2700  Merced  Street 
San  Leandro,  CA  94577 

1  Radkowski  Associates 

ATTN:  Peter  R.  Radkowski 

P.O.  Box  5474 
Riverside,  CA  92517 

4  R&D  Associates 

ATTN:  Jerry  Carpenter 

J.  G.  Lewis 
Technical  Library 
Allan  Kuhl 
P.0.  Box  9695 
Marina  del  Rey,  CA  90291 

1  RCA  Government  Communications 

Systems 

13-5-2  Front  &  Cooper  Streets 
Camden,  NJ  08102 

2  Science  Applications,  Inc. 

ATTN:  Burton  S.  Chambers 

John  Cockayne 
P.O.  Box  1303 
1710  Goodridge  Drive 
McLean,  VA  22102 

1  Science  Applications,  Inc. 
ATTN:  Technical  Library 

P.O.  Box  2351 
La  Jolla,  CA  92038 

1  Systems  Science  and  Software 
ATTN:  C.  E.  Needham 

P.O.  Box  8243 
Albuquerque,  NM  87198 

1  Systems  Science  and  Software 
ATTN:  Technical  Library 
P.O.  Box  1620 
La  Jolla,  CA  92037 


No .  of 

Copies  Organi/..'  i  i  ■  m 

1  TRW  Systems  Group 

ATTN:  Benjamin  Sussholtz 

One  Space  Park 

Redondo  Beach,  CA  92078 

2  Union  Carbide  Corporation 

Holifield  National  laboratory 
ATTN:  Doc  Control  for  Tech  Lib 

Civil  Defense  Research  Proj 
P.O.  Box  X 

Oak  Ridge,  TN  37830 

1  Weidlinger  Assoc.  Consulting 
Engineers 

ATTN:  M.  L.  Baron 

110  East  59th  Street 
New  York,  NY  10022 

1  Battelle  Memorial  Institute 
ATTN:  Technical  Library 

505  King  Avenue 
Columbus,  OH  43201 

1  California  Inst  of  Tech 

ATTN:  T.  J.  Ahrens 

1201  E.  California  Blvd  . 
Pasadena,  CA  91109 

2  University  of  Denver 

Denver  Research  Institute 
ATTN:  Mr.  J.  Wi sot ski 

Technical  Library 
2390  S.  University  Blvd 
Denver,  CO  80210 

1  IIT  Research  Institute 
ATTN:  Milton  R.  Johnson 

10  West  35th  Street 
Chicago,  IL  60616 

1  J.  D.  Haltiwanger 
Consulting  Services 
B106a  Civil  Engineering  Bldg. 
208  N.  Romine  Street 
Urbana,  IL  61801 


15 


DISTRIBUTION  LIST 


No .  of  No .  o f 

Copies  orj-.m  i  z.-it  ion  Copies  Organ  i  zn  t  i  on 


)  Massarhnset l s  Institute  of 
Technology 

Ae roe  list i e  and  Structures 
Research  laboratory 
ATTN:  Dr.  K.  A.  Witmer 

Cambridge,  MA  02119 

P  Southwest  Research  Institute 
ATTN:  Dr.  W.  K .  Baker 

A.  B.  Wenzel 
8500  Culehru  Road 
San  Antonio,  TX  78206 


Dir,  CSAMSAA 

ATTN:  DRXSY-I) 

DKXSY-MP,  H 
Cdr,  USATKCOM 

ATTN:  DRSTF.-TO-F 

Dir,  USACSL 

Bldg.  F3516,  FA 
ATTN:  DRDAR-CI.B 


1  SRI  International 

ATTN:  Dr.  0.  R.  Abrahamson 

313  Ravenswood  Avenue 
Menlo  Park,  CA  94025 


1  Stanford  University 

ATTN:  Dr.  D.  Be r shader 

Durand  Laboratory 
Stanford,  CA  94305 


1  Washington  State  University 
Physics  Department 
ATTN:  G.  R.  Fowles 

Pullman,  YA  99164 


Cohen 


PA 
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USER  EVALUATION  OF  REPORT 


i’lease  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  place 
in  the  mail.  Your  comments  will  provide  us  with  information  for 
improving  future  reports. 

1 .  BRL  Report  Number__ _ 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc . ) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc . ) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


Director 

US  Army  Ballistic  Research  Laboratory 
ATTN:  DRDAR-TSB 

Aberdeen  Proving  Ground,  MD  21005 


FOLD  HERE 


